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SUMMARY 


A method has "been developed for determining the- stahility of a scalar 
neutral equation with constant coefficients and constant time delays. A 
neutral equation is hasically a differential equation in which the highest 
derivative appears hoth with and without a time delay. Time delays may appear 
also in the lower derivatives or the independent variahle itself. The method 
is easily implemented and an illustrative example is presented. 


INTRODUCTION 


Ordinary differential eqiiations with time delays are called differential- 
difference equations (ref. l). Two basic types of differential-difference 
equations are retarded and neutral equations. The stability of the solutions of 
these equations is related to the roots of a characteristic equation* Generally 
this characteristic equation is transcendental and thus has an infinite number 
of roots. 

A convenient method is developed in reference 2 for examining the 
stability of retarded equations with many time delays (not necessarily distinct) 
and a scalar neutral equation with one delay. The purpose of the present paper 
is to develop the' basic method of reference 2 for neutral equations with many 
time delays. 
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Mathematical notations: 


I [ absolute value or magnitude 

arg argument 

0^ arbitrarily small positive values 

Dots over a symbol denote derivatives with respect to time . 

ANALISIS 

A method is developed herein for determining the stability of the neutral 
equation 

1 [a + b x^^^(t - T.)] = 0 (1) 

• 3=0 ^ 

where a^g ^ 0^ bj^ ^ 0^ 0 < ^ for 0“ O^l.^*^* y H* 1^ 

and x^^nt ) denotes the j'th derivative of x(t). 

The characteristic equation associated with equation (l) is 


H 

L(s) = S 
3=0 


(a . + b . e 


-T . 
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( 2 ) 


It has been shown (ref. 3) that if all the roots s = g- + in> of 
aquation (2) satisfy the property 


g ^ ~ iu < 0 


(3) 


(There y is a positive constant^ then the solution of equation (l) is of 
jxponentlal order as t -♦ qo; that is 

]x(t)f < d e”®^ (4) 

There d > 0 is a constant real number and c is arbitrary on the 
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interval (O, l). Hence, if all the characteristic roots have negative real 
parts and are not asymptotic to the imaginary axis, then x(t) -♦ 0 as t ^ os 
( asymptotically stable ) . 

If there is a root of L(s) with positive real part, then eq.uation (l) 
has a divergent mode and is said to be -unstable. 


Helative Stability 

If it can be determined that there are no roots of the characteristic 
equation with real parts greater than a specified negative real number, then 
the solution to the neutral equation is asymptotically stable. 

Relative stability for a specified value u of ju in equation (3) is 
indicated herein by the number of roots of the characteristic equation with 
Q > -/?. For example, the neutral system is said to be relatively more stable 
when all the roots satisfy q- < -ju < 0, than when there is a root with 
-U < q < 0. Relative stability boundaries in the plane of two system parameters 
are boundaries corresponding to a root with q = -U* 

The stability method to be presented is based on determining the number of 
roots of the characteristic equation with real parts greater than a specified 
negative real number The method is convenient for determining the number 

of roots of the characteristic equation with real parts located between specifie 
negative real mmibers. The approach consists of separately examining the 
arbitrarily large modulus roots and the finite roots. The large modulus roots 
are examined by using a simple ex-pression for their asymptote; whereas, the 
finite roots are examined by computing the magnitude of a complex-valued functio 
on a finite interval. 


large Modulus Roots 


All roots of equation (2) must satisfy the inequality 






s!^^ Y | fa.[ + lb|! 
3=0 




(5) 


obtained from equation (2). It can be shown that since a^ ^ 0 and 4 0, 
the roots have bo-unded q. Hence, in order for the large modulus roots 
(Is I -♦ oo) to satisfy equation (5) 
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From eqimtion (6), g beccffiies arbitrarily close to 


CT = - — In 
“ % 


% 


b. 


N 


( 7 ) 


This relation represents the as 3 nnptote of the large modnlus roots and is showi 
graphically in figvire 1. 


For 
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< 1 in figTore 1, q- > 0 j and equation (3) with o = o’ is not 
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in figure 1. Then, g> = o satisfies equation (3 ) with p = w, whenever 
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( 8 ) 


There are then no infinitely large modulus roots with o- > -a in the neutral 
system. It remains to examine the number of finite roots with o > -u. 


Finite Roots 


For Tj -♦ o'*’ , L(s) has 
J 

the polynomial equation 


F roots arbitrarily close to the F roots of 
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S 

3=0 



0 


( 9 ) 


and the remaining roots have arbitrarily large moduli (ref, 2). For 

^F A 

I >1 in equation (?), e, 

00 


Tjj -♦ o'*" and 


b. 


F 


> 1 in equation (?)> o < -a. Therefore, L(s) and 


Lq(s) have the same number of roots with a > -m (initial relative stability). 
Since the complex roots occur in complex conjugate pairs , only roots with non- 
negative imagimry parts (co ^ O) are considered. 


As one of the time delays, say Tt^ is increased in a continuous manner 
with the remaining delays held fixed, me finite roots of L(s) move in some 
continuous manner (ref. 2), generating root locxxs curves in the complex root 
plane (s-plane). 
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Intersection Points s = -A? + iw and 

Corresponding Delays 

A root locus curve must intersect the -Al-line (dashed line) in figijre 2 
in order for the number of roots of L(s) ■with g* > -M to change. These 
intersection points (-in, cd) and the corresponding values of the delay 
which result in these intersection points are discussed in this section. The 
change in the relative stability as a root locus ctirve crosses an Intersection 
point is presented in the next section. 

For a specific time delay ecjuation ( 2 ) can be written as 

* 'V s 

L(s) = Hjj(s) - Jj^(s) e ^ e 0 (lO) 

where 

H^(s) - S a. s*^ + S b. s‘^ e ^ (ll) 

d=0 ^ J=0 ^ 

and 

Jj^(s) = -bj^ s^ (12) 


At an intersection point s - -ju + im, equation (lO) is equivalent to 

/\ 

Ih^(-U, cjl)) 1 == |jj^(-4, oi))l e 

and 


Tr, = - arg 
K CO 


Jjr(-M, 03 ) 


03 ) 


+ 2prt 


where co) = Hp^(-h + io)), co) - jg(-M + ijo)^ and 


- jt < arg. 


03 ) 

Hp^(“ 4 j 03) 


< rt 


(13) 


( 14 ) 


(15) 


It is assumed that co 7^ 0 and co) 7^ 0 . To handle these special cases, 

the approach used in reference 2 -m-j be followed. Only non-negative values of 
the integer p in equation (l 4 ) are of interest because ^ 0 and co > 0. 
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Equation (13) gives the points (-m, oj) where the root locus curves 
intersect the - 4 - line in figure 2., and equation (l4) gives the corresponding 
values of Tg- which result in these intersection points. In general, the 
values of oj at an intersection point must be fomd by an iteration process. 
The values of m which may satisfy equation (l4) are restricted to some finite 
interval (O, where m is an upper bound on to determined from 

equation (5)* Also, a useful bound on the integer p in equation (l4) is 
obtained as 



1 iU Tt^ 

1 . m K 

-t* W WW 

2 2 * 


(16) 


where 



and CO ^ co^. 


Change in ETumber of Roots 
With q > -iu 

Let 11 ( 715 -, denote the number of roots of L(s) with gr > -ju at Xg- 
for fixed x -, j 7 ^ Kj and define the testing function 

t) 






(17) 


Then, the following theorem can be used to determine the change in the number 
of roots of L(s) with q > as Xg- varies. 

Theorem : Let (- 4 , co) be an intersection point with corresponding delay Xg^. 

Let < CO and Ug > co be real numbers for which Wj^(-ju, and Wj^(- 4 , Ug) 

are defined, and such that there are no other intersection points with 
imaginary parts which lie on the interval [a^ Ug]. Then, for e an 
arbitrarily small positive number 
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|W^(-ju, a^) I and oc,g | are greater than 1 or both less than 1. 

This theorem is developed in reference 3 by extending the t-decomposition 
method, as refined by Lee and Hsu (ref. 4). 


The theorem is interpreted as follows* Let (-iE?, ao) be an intersection 
point, where jU is specified and oj is a root of equation (13)« If this is 
the only value of on the interval ^ ^ a^, which satisfies 

equation (13)^ then the change in the relative stability at the intersection 
point is determined by computing lWg(-u‘, ct,^) [ and |Wg;(-iU, o^ )} . For example, 
from condition 1 of the theorem, if | > 1 and cCg)| <^1, 

then the system gains exactly one root with cr > -m; that is, H(tv + -m)= 

-iu) + 1. 


The values of 


K 


at all the intersection, points are ordered by 


increasing magnitude to obtain the change in the relative stability as 

increases to its final desired value Tj^. As each delay is varied, that 
delay becomes x-^ in the theorem. 


Intersection points (-4, 00 ) satisfy equation (13)^ or fW_(-4, ai)| = 1. 

IX 

In choosing and in the theorem, it is expedient to note that 

|¥j^(-jU, co) ] increases as p increases for each value of cu e(0^ oci^]* 


APPLICATION 


The relative stability of the neutral equation 


,0102l<- t(t) + .00704 t(t) + .250 t(t) + .l 63§t(t - Tj^) = 0 (18) 

where ^ is a system gain constant and v > 0 is a constant time delay may 

now be determined. This equation was used in reference 5 in examining a yaw 
damper control system for an airplane with rudder deflection made proportional 
to the yawing acceleration. 

The characteristic equation associated with eqiiation (18) can be written 
as 


■•T S 

L(s) = P(s) - Q(s)fe ^ = 0 

where 

P(s) = .01024s^ + .00704s + .250 
Q(s) = .163 s^ 


(19) 

( 20 ) 
( 21 ) 
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With s = -ju + ia>) equation. (I 9 ) can be used to -write 


and 



(22) 


(23) 


Ifow, with u specified, equations ( 22 ) and ( 23 ) can be used to -partition the 
plane of f and into differen-b regions as oj > 0 is allowed to -vary. 
The solid lines in figure 3 were generated in this manner. Any point on a 
partitioning line or boundary corresponds to a root locus curve intersecting 
the -n-line in figure 2 . 


To examine the stability condition (stable or -unstable) or the number of 
roots with o > -w in the regions of figure 3 # it is useful to -write 
equation (I 9 ) in the form 


•"X s 

L(s) = Hg.(s) - Jj^(s)e ^ =0 (2k) 

where 

= p(s) (25) 

and 


Jjr(s) = §.Q(s) 


( 26 ) 


The initial stability of equation (24) along the ^-axis (xg. -♦ 0"^) is 
e-valuated by using equations ( 7 ) and ( 9 ), which became 


^ ~ In 


,0102k 

. 1 ^ 3 ^ 


and 


( 2 ?) 


Lq(s) = (.01024 + .l63f)s^ + .00?04s + .250 = 0 (28) 

For -» 0’*’ and ^ = .04, there is one root -with g = -.21 and g- 

As Tg- increases frcm O"*" with ^ = .04 in figure 3> the relative°°stability 
boundary for -u = -. 5 is. intersected. For this intersection point, it can 
be shown that and in the theorem can be chosen as a, = 3 and 

ttg = 4. Then, since |¥g.(-.5, 3)|< 1 and |Wg-(-.5, 4)| >1, condition 2 of 
the theorem applies. Thus, the neutral system loses one roo-t with 0 > -. 5. 
(This is the root which originally had g- = -.2L) Inside the closed region for 
-4 = -.5, there are no roots -with g > -. 5. This same procedure is \osed to 
determine which side of the curves in figure 3 should be hatched. 
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The hatching convention is as foUovrs: Ikssiiag from the hatched (nnhatched) 

side of a hountory line ccaresponding to a particular value of M to the 
unhatched (hatched) side of the boundary results in the gain (loss) of exactly 
one root with o > -jE?. 

At the point (§, Tg) = (.04, .2), the system has no roots with o > -.J* 

The value of a in equation ( 2 j) at this point is g* = -2,257, 

00 00 


CaSCEUDIlG REMARKS 


A method has been developed fen* determining the stability ani relative 
stability of scalar neutral equations, with constant coefficients and constant 
time delays. The approach was to determine the number of roots of the 
characteristic equation with real parts greater than specified negative real 
numbers. The method consists of separately examining the large modulus roots 
and finite roots. The large modult^ roots are examined by t^ing a simple 
expression for their asymptote; the finite roots are examined by computing 
the nagnitude of a complex-valued function on a finite Interval. 

The stability method is convenient for determining the number of roots of 
the characteristic equation with real parts located between specified negative 
real numbers. An example which has occurred in practical application has been 
provided to illustrate the method. 
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